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SUMMARY 


The  generation  of  long  period  (20-50  second)  surface  waves  by  an 
explosion  on  an  island,  inside  a  mountain,  or  near  a  material  discontinuity  is 
studied  using  two-dimensional  axisymmetric  finite  difference  calculations.  The 
generation  of  surface  waves  can  be  reduced  substantially  if  the  explosion  is 
close  to  a  vertical  boundary  such  as  a  material  velocity  reduction,  mountain 
slope,  or  island/ocean  boundary.  This  occurs  because  surface  waves  from  a 
shallow  explosion  are  generated  primarily  by  the  horizontal  stress  components 
which  are  reduced  at  the  boundary. 

Three  sets  of  finite  difference  calculations  were  performed  for  an  "island" 
surrounded  by  an  "ocean"  of  air,  water,  and  a  low  velocity  solid.  The  surface 
waves  are  measured  on  "land"  on  the  other  side  of  the  "ocean."  A  "solid* 
calculation  with  a  uniform  structure  is  also  performed  for  comparison. 
Calculations  include  "ocean"  depths  of  3  km  and  6  km,  and  the  "island"  varies  in 
diameter  from  6  km  to  48  km.  For  a  48  km  island,  the  material  boundary  has  little 
effect  on  the  surface  waves.  For  an  island  width  equal  to  the  ocean  depth, 
however,  the  surface  wave  generation  is  reduced  by  an  order  of  magnitude  for 
air  and  water,  and  reduced  by  more  than  a  factor  of  2  for  the  low  velocity  solid. 
Comparisons  with  the  geography  of  the  Amchitka  and  Mururoa  island  test  sites 
show  that  this  effect  is  small  for  explosions  at  those  sites,  however  the  effect 
appears  to  be  important  for  Novaya  Zemlya  and  other  sites  where  explosions 
are  detonated  inside  mountains,  and  at  sites  with  strong  material  variations. 

To  estimate  the  effect  of  a  high  velocity  boundary  surrounding  a  low 
velocity  source  region,  we  performed  a  calculation  for  an  explosion  in  a 
sedimentary  basin  modeled  after  Yucca  Valley  at  the  Nevada  Test  Site.  In  order 
to  include  finer  details  of  the  source  region,  a  finely  gridded  finite  difference 
calculation  was  performed  in  the  source  region,  and  surface  waves  exterior  to 
the  calculation  were  calculated  using  the  representation  theorem.  We  find  that 
surface  waves  are  amplified  by  about  a  factor  of  two  by  the  high  velocity 
interface  surrounding  the  low  velocity  valley. 


1.  INTRODUCTION 


We  investigate  the  effect  of  a  near-source  boundary  or  material  property 
change  on  the  generation  of  long  period  surface  waves  by  explosions.  This 
work  was  initially  motivated  by  the  work  of  von  Seggern  (1978),  and  others  who 
found  a  possible  bias  between  Amchitka  and  NTS  and  suggested  that  surface 
waves  generated  by  an  explosion  on  an  island  might  differ  from  those 
generated  by  a  mid-continent  explosion.  The  surface  wave  amplitudes  from  the 
Amchitka  explosion  Milrow,  for  example,  were  significantly  lower  than  the 
surface  wave  amplitudes  for  the  NTS  explosions  Boxcar,  Benham,  and 
Handley,  all  of  which  had  similar  yields  (von  Seggem,  1978;  Stevens,  1986b). 
Earlier  investigations  have  looked  at  the  effect  of  tectonic  release  (Toksoz  and 
Kehrer,  1972)  and  large  scale  three-dimensional  earth  structure  (McLaughlin,  et 
al,  1991). 

There  is  a  common  perception  that  surface  waves  should  not  be  affected 
strongly  by  the  details  of  the  source  region,  because  surface  waves  are  long 
period  phenomena  with  wavelengths  much  larger  than  a  test  site.  However, 
surface  waves  exist  because  of  the  free  surface  boundary  condition,  and 
changes  in  the  boundary  shape  can  affect  surface  wave  generation.  This  is 
particularly  true  for  explosions,  because  long  period  surface  waves  from 
explosions  are  generated  by  the  static  horizontal  stress  field  that  remains  after 
the  explosion  has  been  detonated.  If  the  stress  field  is  relieved  by  the  presence 
of  a  vertical  boundary  such  as  an  island/ocean  interface  or  the  side  of  a 
mountain,  then  the  long  period  surface  waves  can  be  sharply  reduced,  even 
though  the  scale  of  the  structural  changes  are  much  smaller  than  a  wavelength. 
In  addition  to  Amchitka,  the  French  nuclear  test  site  at  Mururoa  is  located  on  an 
island  that  drops  very  steeply  into  the  South  Pacific  ocean.  At  the  old  French 
Sahara  test  site  and  the  two  Soviet  test  sites  at  Novaya  Zemlya  and  Degelen 
Mountain,  explosions  are  detonated  inside  mountains. 

We  have  performed  a  series  of  finite  difference  calculations  to  model  the 
effect  of  the  island/ocean  interface  on  explosion-generated  surface  waves. 
"Island"  is  used  here  in  a  generic  sense  to  mean  a  source  region  of  higher 
velocities  than  the  surrounding  regions.  We  have  performed  three  sets  of 
calculations:  first  with  the  island  surrounded  by  air;  second  with  the  island 


surrounded  by  water;  and  third  with  the  island  surrounded  by  a  low  velocity 
solid.  The  calculations  therefore  apply  not  only  to  an  explosion  on  an  island,  but 
also  to  an  explosion  in  a  mountain,  or  to  an  explosion  in  a  high  velocity 
intrusion  such  as  a  salt  plug  surrounded  by  low  velocity  sediments.  We  look  in 
particular  at  the  effect  of  the  size  of  the  island  and  the  depth  of  the  ocean  on  the 
surface  wave  amplitudes.  We  show  in  the  following  sections  that  the  surface 
wave  amplitudes  can  be  reduced  significantly  by  the  presence  of  a  vertical 
boundary  near  the  source. 

The  calculations  described  above  were  performed  with  a  fairly  coarse 
grid  spacing  in  order  to  directly  calculate  surface  waves  several  hundred 
kilometers  from  the  explosion.  In  order  to  better  quantify  the  effect  of  depth  on 
surface  wave  generation,  we  also  performed  a  series  of  plane  strain 
calculations  of  the  dilatation  at  several  depths  due  to  a  Rayleigh  wave  incident 
on  a  mountain  and  then  used  reciprocity  to  determine  the  effect  of  depth  on 
Rayleigh  wave  generation  by  explosions  at  these  depths.  We  find  that  the 
Rayleigh  wave  amplitude  is  reduced  significantly,  even  for  an  explosion  at  the 
base  of  the  mountain,  and  is  strongly  reduced  for  explosions  within  the 
mountain. 

We  also  performed  an  axisymmetric  calculation  of  an  explosion  in  a 
mountain  modeled  after  a  profile  of  a  mountain  at  the  Novaya  Zemlya  test  site. 
In  this  case,  we  used  a  much  more  finely  gridded  model  than  in  the  calculations 
discussed  above  and  used  the  representation  theorem  to  generate  surface 
waves  at  distances  outside  the  grid.  This  was  done  by  saving  the  displacements 
and  stresses  on  a  monitoring  surface  in  the  exterior  structure  and  then 
integrating  together  with  a  multi-modal  surface  wave  Green's  function  to 
generate  surface  waves.  Explosions  were  located  within  the  mountain  and  just 
below  the  base  of  the  mountain.  The  results  agree  with  the  plane  strain 
reciprocal  experiment  discussed  above.  The  surface  wave  amplitudes  are 
reduced,  and  reduced  substantially  more  for  the  explosion  within  the  mountain. 

We  investigate  the  effect  of  a  high  velocity  boundary  on  the  generation  of 
surface  waves  by  an  explosion  in  a  low  velocity  medium.  A  two-dimensional 
model  was  developed  for  Yucca  Valley  at  the  Nevada  Test  Site,  which  is  a 
region  of  low  velocity  sediments  surrounded  by  higher  velocity  rock.  The 


calculation  used  a  finely  gridded  source  region  and  used  the  representation 
theorem  to  generate  surface  waves,  as  with  the  Novaya  Zemlya  mountain 
calculation.  The  resulting  surface  waves  were  compared  with  synthetic 
seismograms  generated  by  an  explosion  in  a  plane-layered  structure 
corresponding  to  the  structure  at  the  source  and  then  propagating  into  the 
external  plane-layered  structure,  using  a  transmission  coefficient  based  on 
conservation  of  energy.  In  this  case,  we  find  that  the  Rayleigh  wave  amplitudes 
are  substantially  enhanced  by  the  high  velocity  boundary. 


2. 


ISLAND/OCEAN  FINITE  DIFFERENCE  CALCULATIONS 


All  calculations  were  performed  with  the  two-dimensional  axisymmetric 
finite  difference  code  TRES-2D  running  on  a  CRAY  2  supercomputer.  The 
explosion  was  modeled  as  an  isotropic  moment  tensor  source  with  a  time 
function  proportional  to  a2te*at,  a  =  0.1  Hz  (corresponding  to  a  minimum  phase 
source  with  corner  frequency  at  0.1  Hz).  The  source  was  converted  to  a 
pressure  applied  to  two  vertical  zones  on  the  axis  of  symmetry.  The  source  is 
therefore  shaped  as  a  cylinder  with  equal  height  and  radius,  and  is 
approximately  isotropic.  For  the  series  of  island/ocean  calculations,  we  used  a 
uniform  grid  spacing  of  1 .5  km,  which  is  designed  to  allow  accurate  resolution  of 
surface  waves  in  the  20-50  second  period  band.  The  source  was  located 
immediately  below  the  free  surface.  A  schematic  diagram  of  the  geometry  used 
in  the  calculations  is  shown  in  Figure  1 .  Calculations  were  performed  for  island 
radii  ranging  from  3  km  to  24  km  and  for  ocean  depths  of  3  km  and  6  km. 

The  ‘ocean*  in  this  model  is  a  region  with  material  properties  different 
from  the  background  model.  The  background  model  is  a  plane-layered  model 
(listed  in  Table  1)  and  is  equivalent  to  a  north-south  cross-section  near  the 
Amchitka  test  site  in  the  three  dimensional  model  of  McLaughlin,  et  al.  (1991). 
One  calculation  was  performed  as  a  reference  using  the  background  structure 
without  the  ocean.  Three  sets  of  calculations  were  performed  with  different 
material  models  in  the  "ocean."  In  the  f:rst  set  of  calculations,  the  ocean  was 
filled  with  water  using  standard  water  density  and  compressional  velocity.  In  the 
second  set  of  calculations,  a  free  surface  boundary  condition  was  used  on  the 
bottom  and  sides  of  the  ocean,  thus  simulating  air  or  vacuum  in  the  "ocean."  In 
the  third  set  of  calculations,  the  ocean  was  filled  with  a  solid  with  a 
compressional  velocity  of  2.4  km/sec,  a  shear  velocity  of  1.3  m/sec,  and  a 
density  of  1.9  g/cm3.  These  values  are  typical  of  the  volcanic  tuffs  found  at  the 
Nevada  Test  Site.  This  model  therefore  simulates  a  high  velocity  intrusion  into  a 
lower  velocity  medium 


Table  1 

Background  Structure  for  the  Axisymmetrlc  Calculations 


Depth 

km 

Thickness 

km 

P-Velocity 

km/sec 

S-Velocity 

km/sec 

Density 

gm/cc 

6.0 

6.0 

6.90 

3.90 

2.29 

18.0 

12.0 

6.90 

3.90 

2.29 

54.0 

36.0 

8.00 

4.50 

3.30 

78.0 

24.0 

8.00 

4.45 

3.30 

126.0 

48.0 

8.00 

4.40 

3.25 

150.0 

24.0 

8.04 

4.40 

3.37 

186.0 

36.0 

8.02 

4.40 

3.37 

222.0 

36.0 

8.01 

4.40 

3.36 

2C8.0 

36.0 

8.60 

4.66 

3.45 

294.0 

36.0 

8.65 

4.68 

3.46 

330.0 

36.0 

8.73 

4.70 

3.49 

366.0 

36.0 

8.79 

4.73 

3.51 

402.0 

36.0 

8.87 

4.76 

3.53 

o© 

oo 

9.15 

4.95 

3.73 

In  Figure  2,  we  show  the  seismograms  calculated  for  the  uniform 
background  model  together  with  seismograms  for  the  1 .5  km  and  3  km  source 
in  a  6  km  radius/6  km  high  island  surrounded  by  air  at  a  distance  of  390  km  from 
the  source.  The  seismograms  differ  substantially  in  amplitude,  with  the  island 
seismograms  lower  in  amplitude  than  the  background  seismogram. 

In  order  to  examine  the  frequency  dependence  of  the  seismograms,  the 
amplitude  and  phase  of  the  surface  waves  at  the  receiver  were  recovered  by 
narrow  band  filtering  and  phase-matched  filtering  as  described  by  Stevens 
(1986a).  The  source  function  was  deconvolved  during  the  processing,  so  the 
spectra  correspond  to  a  source  with  a  step  function  in  moment.  In  Figure  3,  we 
show  the  spectra  between  .01  and  0.1  Hz  for  the  uniform  background  structure, 
and  for  the  water,  air,  and  solid  filled  ocean  cases  with  an  island  radius  of  24 
km,  and  an  ocean  depth  of  3  km.  For  this  small  aspect  ratio  (aspect  ratio  is 
defined  here  to  be  the  ratio  of  the  island  height  to  the  island  radius),  the  island 
structure  has  little  effect  on  the  surface  wave  amplitudes.  The  material  property 
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Seismograms  at  a  range  of  390  km  for  the  3  km  radius  and  24  km 
radius,  6  km  high  isiand.  The  surface  wave  is  reduced  dramatically 
for  the  source  in  the  small  island. 
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Spectra  of  surface  waves  recorded  at  390  km  from  a  24  km  radius 
island  across  the  uniform  background  structure,  and  across  an 
"ocean"  filled  with  water,  air,  and  a  lower  velocity  solid.  In  each 
case,  the  "ocean"  is  3  km  deep.  The  "ocean"  material  has  some 
effect  on  shorter  period  surface  waves,  but  there  is  little  difference 
between  the  three  cases  at  periods  greater  than  20  seconds. 


differences  cause  some  variation  in  spectral  amplitude  at  higher  frequencies, 
but  at  frequencies  less  than  0.05  Hz,  there  is  very  little  difference  between  the 
surface  wave  amplitudes  for  any  of  the  four  structures. 

In  Figure  4,  we  show  the  surface  wave  amplitudes  for  an  island  radius  of 
6  km  and  an  ocean  depth  of  6  km.  In  this  case  the  amplitudes  are  very  strongly 
reduced  for  all  three  ocean  types.  At  25  second  period  (0.04  Hz),  the  amplitudes 
are  reduced  relative  to  the  uniform  case  by  factors  of  2,  6,  and  20  for  the  solid, 
water,  and  air  filled  ocean  cases,  respectively.  In  Figures  5,  6,  and  7,  we  show 
the  25  second  amplitude  as  a  function  of  the  aspect  ratio  of  the  island.  In  each 
case,  the  amplitude  is  reduced  for  aspect  ratios  greater  than  about  0.2,  and  the 
amplitudes  are  reduced  by  a  factor  of  2  at  aspect  ratios  of  0.35,  0.4,  and  0.7  for 
air,  water,  and  solid,  respectively.  For  the  solid  model  used  here,  the  amplitude 
is  reduced  by  a  maximum  of  about  a  factor  of  2.5,  but  the  amplitude  reduction 
can  be  much  larger  for  the  air  and  water  cases. 

In  Figure  8,  we  show  the  bathymetry  near  the  Amchitka  and  Mururoa  test 
sites.  The  Amchitka  bathymetry  was  derived  by  drawing  a  line  southwest  from 
the  test  site  using  maps  from  Olsen,  et  al.  (1972).  The  Mururoa  bathymetry  was 
derived  similarly  by  drawing  a  line  to  the  northwest  of  the  test  site  using  maps 
from  DMA  (1983).  For  Amchitka,  the  aspect  ratio  is  only  about  0.05,  so  clearly 
the  amplitude  reduction  from  the  island  structure  at  Amchitka  is  negligible.  At 
Mururoa,  however,  the  ocean  drops  off  much  more  quickly  away  from  the  test 
site,  and  the  aspect  ratio  of  the  island  is  approximately  0.25.  So  surface  waves 
might  be  reduced  slightly  at  the  Mururoa  test  site  by  this  effect. 
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Figure  4.  Spectra  of  surface  waves  recorded  at  390  km  from  a  6  km  radius 
island  across  a  6  km  deep  "ocean"  of  water,  air,  and  solid.  The 
surface  wave  amplitudes  are  sharply  reduced  at  all  frequencies  in 
all  cases. 
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Figure  6. 


25  second  surface  wave  amplitude  as  a  function  of  the  aspect  ratio 
for  an  air-filled  "ocean." 
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Figure  8.  Bathymetry  near  the  Amchitka  and  Mururoa  test  sites. 
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3.  SURFACE  WAVES  FROM  AN  EXPLOSION  IN  A  MOUNTAIN 


Several  nuclear  test  sites  are  in  areas  characterized  by  mountains  and 
steep  slopes.  Explosions  at  Novaya  Zemlya,  Degelen  Mountain,  and  the  old 
French  Sahara  test  site  have  all  been  located  within  and  close  to  mountains. 
The  results  of  the  previous  section  suggest  that  this  topography  could  have  a 
strong  effect  on  surface  wave  generation.  An  important  issue  is  the  effect  of 
depth  within  and  below  a  mountain  on  surface  wave  generation.  That  is,  will  the 
predicted  amplitude  reduction  vanish  if  the  explosion  is  below  the  base  of  the 
mountain,  or  will  the  distortion  of  the  local  stress  field  by  the  non-planar  surface 
continue  to  cause  an  amplitude  reduction?  The  numerical  tests  described  in  the 
previous  section  used  a  grid  too  coarse  to  answer  this  question. 

In  order  to  examine  this  effect  in  more  detail,  we  have  performed  two 
numerical  experiments.  First,  we  performed  a  series  of  plane  strain  calculations 
for  a  Rayleigh  wave  incident  on  a  mountain,  using  reciprocity  to  determine  the 
Rayleigh  wave  generated  by  an  explosion  in  the  mountain  from  the  dilatation 
produced  by  the  Rayleigh  wave.  Second,  we  performed  a  calculation  for  an 
explosion  in  a  mountain  modeled  after  the  actual  topography  of  a  mountain  at 
the  Novaya  Zemlya  test  site. 

3.1.  Reciprocity  Experiment 

The  geometry  of  the  plane-strain  reciprocal  experiment  is  shown  in 
Figure  9.  A  planar  Rayleigh  wave  pulse  is  incident  from  the  left  upon  a  steep 
mountain  with  45  degree  slopes  and  height  L=  12  units  and  radius  of  16 
(aspect  ratio  0.75).  The  half-space  had  a  compressional  velocity  of  6  units/sec, 
and  a  shear  velocity  of  3.55  units/sec.  Vertical  displacements  and  dilatations 
are  recorded  at  locations  indicated  on  the  figure.  Incident  displacement  and 
dilatation  are  then  used  to  normalize  each  synthetic  seismogram  to  the  incident 
pulse.The  seismogram  locations  correspond  to  the  top,  middle,  and  bottom  of 
the  mountain  as  well  as  a  location  deep  under  the  mountain  (12  units)  and  a 
location  on  the  far  side  of  the  mountain  (transmitted)  at  two  depths  (deep  and 
shallow). 
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Geometry  of  the  2-D  plane  strain  reciprocal  experiment.  The 
vertical  displacements  and  dilatations  due  to  an  incident  Rayleigh 
wave  were  measured  at  the  locations  shown.  This  is  equivalent  to 
measuring  the  far  field  Rayleigh  wave  generated  by  a  vertical 
point  force  and  explosion,  respectively. 


By  reciprocity  the  ratio  of  the  recorded  displacement  to  the  incident 
Rayleigh  wave  vertical  displacement  is  the  equivalent  to  the  far-field  Rayleigh 
wave  displacement  from  a  vertical  point  source.  Similarly,  the  ratio  of  the 
recorded  dilatation  to  the  incident  Rayleigh  wave  is  equivalent  to  the  far-field 
Rayleigh  wave  from  a  dilatational  (explosion)  source.  Therefore,  we  can 
examine  the  relative  far-field  Rayleigh  wave  excitation  from  sources  within  and 
under  the  mountain.  The  effect  of  transmission  through  the  mountain  from 
source  to  receiver  can  be  examined  as  well. 


Figure  10  shows  the  spectral  amplitude  ratio  of  the  vertical 

displacements  to  the  incident  Rayleigh  wave  displacement  as  a  function  of 

l 

normalized  wavenumber,  k' =  — — ,  where  L=  12  units  is  the  height  of  the 

V 

mountain,  c  =  3.252  units/sec  is  the  Rayleigh  wave  phase  velocity,  and  f  is 
frequency.  This  is  the  reciprocal  problem  of  the  Rayleigh  wave  generated  by  a 
point  force  at  the  locations  shown.  We  consider  k'  «  1  to  be  the  long  period 
limit  while  k'  -  1  corresponds  to  the  resonances  in  the  scattering  problem. 
Note  that  the  effect  of  a  vertical  point  force  at  the  top,  middle,  bottom,  or  deep 
under  the  mountain  is  not  much  different  than  a  source  to  the  side  of  the 
mountain  with  transmission  through  the  mountain.  Some  amplitude  is  lost  due 
to  scattering  but  from  k'=0.1  to  k'=1,  the  effect  is  generally  less  than  30  percent. 
For  k'<1,  the  effect  of  the  mountain  for  a  vertical  point  force  is  comparable  to  the 
transmission  through  the  mountain  and  differs  by  less  than  30  percent. 

However,  if  we  examine  the  problem  of  explosive  sources  at  these 
locations,  the  results  are  quite  different.  Figure  11  shows  spectral  amplitude 
ratios  for  dilatational  recordings  at  the  respective  source  locations  normalized  to 
the  incident  Rayleigh  wave  dilatation.  Note  that  a  source  deep  under  the 
mountain  is  much  the  same  as  a  source  to  the  side  of  the  mountain  with 
transmission  through  the  mountain,  and  for  k'=0.1  the  effect  is  less  than  30 
percent.  However,  the  Rayleigh  wave  excitation  from  a  source  at  the  top  of  the 
mountain  is  reduced  by  nearly  90  percent  at  low  frequencies.  The  deeper  the 
source,  under  the  mountain,  the  less  the  effect,  but  a  source  located  at  the 
bottom  of  the  mountain  is  still  reduced  by  nearly  50  percent  with  respect  to  the 
half-space  source  at  low  frequencies. 
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Figure  10.  Spectral  amplitude  ratio  of  the  vertical  displacements  to  the 
incident  Rayleigh  wave  displacement  as  a  function  of  normalized 
wavenumber.  This  demonstrates  that  the  Rayleigh  wave 
generated  by  a  vertical  point  force  is  insensitive  to  the  location  of 
the  point  force. 
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Figure  11.  Spectral  amplitude  ratios  for  dilatational  recordings  at  the 
respective  source  locations  normalized  to  the  incident  Rayleigh 
wave  dilatation.  This  demonstrates  that  the  Rayleigh  wave 
amplitude  generated  by  an  explosion  is  quite  sensitive  to  the 
location  of  the  explosion  within  the  mountain.  The  Rayleigh  wave 
amplitude  is  sharply  reduced  by  location  in  or  immediately  under  a 
mountain. 


This  effect  has  significant  implications  for  long-period  Rayleigh  wave 
excitation  from  underground  explosions  placed  within  a  mountain.  The 
perturbation  of  the  free  surface  alters  the  long-period  response  of  the  medium  to 
a  dilatational  source  much  more  than  a  vertical  force.  The  reciprocal  statement 
is  that  the  free  surface  perturbation  alters  the  strain  field  response  much  more 
than  the  displacement  response  for  an  incoming  wavefield. 

3.2.  Novaya  Zemlya  Mountain  Calculation 

The  rugged  topography  of  the  Novaya  Zemlya  test  site  is  shown  in  a 
bedsheet  plot  based  on  Defense  Mapping  Agency  (DMA)  data  in  Figure  12. 
Mountains  rise  from  sea  level  to  over  one  kilometer  in  the  region  where  most 
tests  have  been  located  (based  on  locations  from  Lilwall  and  Marshall,  1986). 
Higher  elevations  and  steeper  topography  occur  outside  the  region  where  tests 
have  occurred  to  date.  In  Figure  13,  we  show  a  SPOT  satellite  image  displayed 
with  the  DARPA  Yield  Estimation  System  (Murphy,  et  al,  1991)  showing  the 
central  part  of  the  Novaya  Zemlya  test  site  near  the  Matochkin  Shar  strait 
showing  the  detailed  local  topography  in  this  ro"«on. 

A  profile  was  constructed  using  the  actual  topography  of  one  of  the 
mountains  in  the  central  Novaya  Zemlya  area.  We  selected  one  of  the  steepest 
slopes  in  order  to  get  an  estimate  of  the  maximum  effect  that  such  topography 
could  have  on  surface  wave  generation.  This  was  a  profile  running  to  the 
southeast  from  the  highest  peak  in  the  central  Matochkin  Shar  region.  The 
profile  was  then  turned  into  a  finite  difference  grid  with  grid  elements  100 
meters  by  100  meters.  The  surface  profile  of  the  finite  difference  grid  is  shown  in 
Figure  14. 

The  grid  spacing  required  by  this  calculation  was  too  fine  to  extend  to  the 
distances  of  interest  for  measuring  surface  waves.  Instead,  we  saved  the 
displacements  and  stresses  on  a  monitoring  surface  outside  the  source  region 
(below  and  outside  of  the  elevated  region),  and  used  the  representation 
theorem  to  calculate  seismograms  at  the  desired  distances.  This  technique  has 
been  described  in  detail  by  Bache,  et  al.  (1982),  Day,  et  al.  (1983),  and 
Stevens,  et  al.  (1991).  The  finite  difference  calculation  was  performed  with 
uniform  material  properties  throughout  the  grid  (P  velocity  6.0  km/sec,  S  velocity 
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Figure  12.  Bedsheet  plot  showing  the  topography  of  the  Soviet  Novaya 
Zemlya  test  site.  The  rectangle  indicates  the  region  where  most 
tests  have  taken  place. 
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Novaya  Zemlya  Mountain  Axisymmetric  Model 
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Figure  14.  Mountain  topography  model  and  explosion  source  locations  used  in  the  finite  difference  calculations  of  the 
Novaya  Zemlya  mountain. 


3.55  km/sec,  density  2.50  gm/cm3).  For  calculation  of  the  surface  wave 
seismograms,  we  used  the  same  structure  used  in  the  previous  set  of 
calculations  (see  Table  1),  but  with  the  upper  two  kilometers  modified  to  be  the 
same  as  those  used  in  the  finite  difference  calculation.  The  source  function 
used  in  the  calculation  was  the  same  as  in  the  previous  calculations. 
Calculations  were  performed  for  depths  of  600  meters  and  1300  meters  below 
the  top  of  the  mountain  as  shown  in  Figure  14.  Since  the  top  of  the  mountain  is 
1000  meters  above  sea  level,  these  depths  correspond  to  locations  inside  the 
mountain  and  just  below  the  mountain.  Calculations  were  performed  for  15 
modes,  although  the  fundamental  mode  Rayleigh  wave  dominates  in  the 
frequency  range  of  interest. 

Synthetic  surface  wave  seismograms  were  generated  from  these  two 
finite  difference  calculations  at  a  distance  of  390  km.  A  synthetic  seismogram 
was  also  calculated  for  a  point  explosion  source  in  the  plane  layered  model 
without  the  mountain  at  a  depth  of  600  meters.  Unlike  the  mountain  simulations, 
the  results  of  the  plane-layered  calculation  are  insensitive  to  depth  for  the  low 
frequencies  of  interest  here.  The  seismograms  were  phase-match  filtered  as  in 
the  earlier  calculations  to  recover  the  surface  wave  spectra.  The  resulting 
spectra  are  shown  in  Figure  15.  The  results  show  that  the  surface  wave 
amplitude  is  reduced  by  almost  a  factor  of  three  for  the  explosion  within  the 
mountain,  and  is  reduced  by  a  factor  of  approximately  1.5  for  the  explosion 
located  just  under  the  mountain  compared  to  the  plane-layered  surface  wave. 
These  results  suggest  that  yield  estimates  based  on  surface  wave  amplitudes 
for  explosions  in  mountainous  regions  such  as  Novaya  Zemlya  could 
underestimate  the  actual  explosion  yield. 


4. 


CALCULATION  OF  AN  EXPLOSION 
IN  A  SEDIMENTARY  VALLEY 


The  previous  sections  were  concerned  with  the  effect  of  a  near-source 
boundary  or  decrease  in  velocity  on  surface  wave  amplitudes.  In  this  section, 
we  examine  the  converse  of  this  problem:  the  effect  on  surface  wave  amplitudes 
of  a  strong  increase  in  velocity  in  the  near-source  region. 

A  2D  axisymmetric  linear  finite  difference  calculation  was  performed  to 
examine  the  excitation  of  surface  waves  from  an  explosion  in  a  deep  low- 
velocity  graben  such  as  Yucca  Flats  at  NTS.  As  with  the  Novaya  Zemlya 
calculation,  the  displacements  and  tractions  were  saved  on  a  monitoring 
surface  outside  the  source  region  and  modal  summation  methods  were  used  to 
propagate  the  seismic  waves  to  a  distance  of  390  km.  This  hybrid  seismogram 
is  referred  to  as  the  representation  theorem  seismogram.  The  finite  difference 
calculation  was  performed  with  a  grid  spacing  of  125  m  and  using  standard 
criteria  of  10  grids  per  wavelength  the  calculation  should  be  good  up  to  0.5  Hz. 
The  source  time  function  was  the  same  function  used  in  the  previous 
calculations. 

The  velocity  model  is  shown  in  Figure  16.  The  axisymmetric  graben  has 
a  1 0  km  diameter  and  is  filled  with  a  250  m  thick  layer  of  low  velocity  alluvium  (a 
=  1 .34  km/s,  3  =  0.603  km/s,  p  =  1 .70  g/cm3  where  a  represents  P  velocity,  p 
shear  velocity,  and  p  density),  a  250  m  thick  layer  of  unsaturated  tuff 
(a  =  2.14  km/s,  p  =  0.963  km/s,  p  =  1.90  g/cm3),  and  a  500  m  thick  layer  of 
saturated  tuff  (a  =  3.00  km/s,  P  =  1 .50  km/s,  p  =  2.00  g/cm3),  The  basement  is 
paleozoic  carbonates  (a  =  5.00  km/s,  p  =  3.00  km/s,  p  =  2.70  g/cm3),  with  a 
surface  weathered  layer  500  m  thick  outside  of  the  graben  (a  =  4.50  km/s,  p  = 
2.25  km/s,  p  =  2.50  g/cm3).  This  model  is  based  on  models  used  by  McLaughlin, 
etal.  (1987)  and  Ferguson  (1981). 

For  the  purposes  of  propagating  waves  outside  the  graben  to  a  distance 
of  390  km,  model  BR3  (Table  2)  was  used.  This  is  a  modified  version  of  the 
basin  and  range  model  BR2  of  McLaughlin,  etal.  (1988).  Plane  layered 
synthetics  were  computed  for  comparison  with  the  representation  theorem 
seismogram.  The  plane  layered  synthetics  were  propagated  in  the  BR3  model 


26 


27 


Figure  16.  Model  used  for  axisymmetric  Yucca  Flats  simulation. 


Table  2 


Depth 

km 

Basin 

Thickness 

km 

and  Range  Model  BR3 

P-Velocity  S-Velocity 
km/sec  km/sec 

Density 

gm/cc 

Q 

0.5 

0.5 

4.50 

2.25 

2.50 

50. 

2.0 

1.5 

5.00 

3.00 

2.70 

3.5 

1.5 

5.80 

3.40 

2.80 

125. 

6.5 

3.0 

6.00 

3.46 

2.90 

125. 

20.5 

14.0 

6.30 

3.64 

2.90 

150. 

24.5 

4.0 

6.50 

3.75 

2.90 

150. 

28.5 

4.0 

6.80 

3.93 

3.00 

150. 

34.5 

6.0 

3.98 

3.00 

150. 

38.5 

4.0 

7.80 

4.50 

3.20 

200. 

42.5 

4.0 

7.84 

4.50 

3.20 

200. 

46.5 

4.0 

7.88 

4.50 

3.20 

200. 

50.5 

4.0 

7.92 

4.50 

3.20 

200. 

54.5 

4.0 

7.96 

4.50 

3.20 

200. 

56.5 

2.0 

4.50 

3.20 

200. 

60.5 

4.0 

8.05 

4.50 

3.20 

200. 

70.5 

10.0 

7.90 

4.40 

3.20 

75. 

80.5 

10.0 

7.80 

4.20 

3.20 

75. 

100.5 

20.0 

7.80 

4.20 

3.20 

75. 

120.5 

20.0 

7.90 

4.50 

3.20 

75. 

180.5 

60.0 

4.55 

3.30 

75. 

oo 

OO 

8.10 

3.40 

150. 

for  explosion  sources  in  the  BR3  model  and  in  the  Yucca  model,  where  the 

•  Yucca  model  was  identical  to  the  BR3  model  below  2  km  depth,  with  properties 
given  by  a  vertical  line  through  the  explosion  source  in  the  Yucca  Valley  model 
in  the  upper  2  km  (Table  3).  In  each  case  the  source  was  a  point  explosion  with 
the  same  reduced  displacement  potential  (\|/«,  not  moment)  as  the 

•  representation  theorem  seismogram.  Propagation  from  the  source  structure  to 
the  path  structure  was  accomplished  using  a  transmission  coefficient  based  on 
conservation  of  energy  (Bache,  eta/,  1978;  Stevens,  1986a).  The  synthetic 
seismograms  were  convolved  with  the  same  source  function  used  in  the 

•  representation  theorem  seismogram.  Twenty-seven  modes  were  computed, 
and  used  in  both  the  representation  theorem  and  the  plane-layered 
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Table  3 

Yucca  Valley  Model 


Depth 

km 

Thickness 

km 

P-Velocity 

km/sec 

S-Veiocity 

km/sec 

Density 

gm/cc 

0.25 

0.25 

1.340 

0.603 

1.70 

0.50 

0.25 

2.140 

0.963 

1.90 

0.50 

3.000 

1.500 

2.00 

1.00 

5.000 

3.000 

2.70 

seismograms,  although  only  the  first  few  modes  have  significant  amplitude  in 
the  seismograms,  and  the  fundamental  mode  dominates  at  frequencies  less 
than  0.1  Hz. 

The  seismograms  are  shown  in  Figure  17.  There  is  a  significant  increase 
in  amplitude  of  the  representation  theorem  seismogram  relative  to  the 
Yucca/BR3  seismogram.  The  spectra  of  the  three  seismograms  are  shown  in 
Figure  18.  The  spectrum  of  the  representation  theorem  seismogram  is 
approximately  twice  the  amplitude  of  the  Yucca/BR3  seismogram.  In  fact,  the 
amplitude  is  intermediate  between  the  spectrum  of  the  Yucca/BR3  seismogram 
and  the  BR3/BR3  seismogram,  in  spite  of  the  fact  that  the  BR3/BR3  seismogram 
had  a  moment  four  times  as  large  as  in  the  representation  theorem  calculation. 

Previously,  Vassiliou,  at  al.  (1987)  investigated  Rayleigh  wave 
propagation  in  a  3D  ultrasonic  scale  model  of  the  Yucca  Flats  graben. 
Ultrasonic  scale  modeling  has  the  advantage  that  3D  models  can  be 
constructed  and  experiments  performed  that  would  be  computationally 
impractical.  They  concluded  that  there  was  no  significant  long  period 
enhancement  of  surface  waves  from  sources  in  a  graben  structure  filled  with 
low  velocity  material.  They  used  a  vertical  point  force  as  a  source  in  the  high 
velocity  medium  outside  the  graben  structure  and  measured  vertical 
displacement  at  locations  in  the  Yucca  Flats  model.  However,  this  was  an 
incorrect  use  of  the  elastodynamic  reciprocity  theorem,  which  requires  the 
experimenters  to  measure  dilational  strain  at  locations  in  the  Yucca  Flats  model 
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Figure  17.  Seismogram  from  the  Yucca  Flats  simulation  (top)  together  with 
plane-layered  seismograms  for  a  source  with  the  same  RDP  in  the 
BR3  model  (bottom),  and  in  the  Yucca  model  transmitted  into  the 
BR3  model  (middle). 
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Figure  1 8.  Spectra  from  the  Yucca  Rats  simulation  and  the  two  plane-iayered 
synthetic  seismograms. 
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rather  than  vertical  displacement.  Consequently,  they  could  only  conclude  that 
vertical  point  force  sources  would  not  exhibit  long  period  enhancement  of 
surface  wave  excitation.  Due  to  the  difficulty  of  measuring  dilatational  strain 
and/or  using  explosion  sources  in  the  physical  scale  models,  it  is  not  possible  to 
truly  model  the  surface  wave  excitation  from  explosions  in  the  3D  model. 

Our  computational  modeling  has  shown  that  strain  is  more  sensitive  to 
lateral  heterogeneity  and  free  surface  topography  than  displacement. 
Consequently  seismic  excitation  from  a  moment  tensor  source  is  more  sensitive 
to  near-source  heterogeneity  than  a  source  composed  of  simple  point  forces. 


5.  CONCLUSIONS 


We  have  completed  a  series  of  finite  difference  simulations  designed  to 
model  the  effect  of  an  island/ocean  interface,  mountain/air  boundary  or  near¬ 
source  material  property  change  on  surface  wave  amplitudes.  The  results 
strongly  suggest  that  explosion  yield  estimates  based  on  surface  waves  may  be 
substantially  affected  by  variations  in  near  source  material  properties,  and  in 
particular  by  near  source  boundaries.  The  effect  of  a  free  surface  or 
ocean/island  interface  close  to  the  source  can  substantially  reduce  surface 
wave  amplitudes.  This  effect  could  be  substantial  for  an  explosion  in  a  salt 
dome  surrounded  by  sediments,  for  example.  A  reduction  in  material  properties 
outside,  but  close  to,  the  source  region  can  also  reduce  surface  wave 
amplitudes  to  a  lesser  extent.  Furthermore,  an  increase  in  material  properties 
near  the  source  region  can  cause  amplification  of  the  explosion-generated 
surface  waves. 

These  results  are  clearly  important  for  seismic  yield  estimation,  and  may 
also  explain  some  anomalies  that  have  been  observed  at  certain  test  sites. 
Some  Novaya  Zemlya  explosions,  for  example,  are  reduced  in  amplitude  at  low 
frequencies  and  phase  shifted  (Stevens  and  McLaughlin,  1988;  Tucker,  et  al., 
1989).  Similarly,  it  has  been  suggested  that  surface  wave  amplitudes  are  higher 
(for  the  same  yield)  for  explosions  in  the  southern  part  of  Yucca  Valley  than  for 
explosions  in  the  northern  part.  Finally,  although  we  have  concentrated  on  the 
effect  of  material  boundaries  on  long  period  surface  waves,  similar  effects  may 
exist  for  regional  surface  waves,  particularly  for  Rg  generated  by  quarry  blasts. 

In  this  study,  we  have  concentrated  on  problems  with  simple  or 
representative  geometries  in  order  to  try  to  quantify  the  general  effects  to  be 
expected  from  the  presence  of  material  boundaries  near  explosions  on  long 
period  surface  waves  under  typical  conditions.  The  same  procedure  could  be 
used,  however,  to  explicitly  model  specific  explosions  or  regions  of  interest  for 
problem  events  that  occur  with  complex  geometries.  The  same  procedures  can 
also  be  applied  to  calculations  of  regional  surface  waves  at  higher  frequencies. 
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